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It is usually assumed that WIMPs interact through spin-independent and spin-dependent 
interactions. Interactions which carry additional powers of the momentum transfer, g^, are assumed 
to be too small to be relevant. In theories with new particles at the ~ GeV scale, however, these 
g'^-dependent interactions can be large, and, in some cases dominate over the standard interactions. 
This leads to new phenomenology in direct detection experiments. Recoil spectra peak at non- 
zero energies, and the relative strengths of different experiments can be significantly altered. We 
present a simple parameterization for models of this type which captures much of the interesting 
phenomenology and allows a comparison between experiments. As an application, we find that dark 
matter with momentum dependent interactions coupling to the spin of the proton can reconcile the 
DAMA annual modulation result with other experiments. 



INTRODUCTION 

The presence of dark matter (DM) in our universe 
is now well established by a variety of astrophysical 
measurements, over a wide range of scales from sub- 
kpc to Gpc. Its clustering and low interaction cross 
sections are supported by the success of the CDM 
framework and the ability of N-body simulations 
to reproduce observed structures. In spite of these 
great successes, we remain ignorant to its detailed 
nature. A direct detection of DM via its recoils 
off of nuclei would confirm its particle nature and 
yield insight into its origin. An examination of the 
recoil spectrum would provide important informa- 
tion about its properties, and possibly the formation 
history of the galaxy. 

If the dark matter is a Majorana fermion x ~ ^ 
supersymmetric neutralino being the most promi- 
nent example - the types of interactions available 
are significantly limited. The dominant scatterings 
are mediated via the operators: 

Osi = {xx)im), (1) 
OsD = {xini5x){qYih(i), (2) 

which give respectively spin-independent (SI) and 
spin-dependent (SD) scattering. As these operators 
typically dominate the interaction rate of Weakly 
Interacting Massive Particles (WIMPs) with nuclear 
targets, direct detection experiments quote results as 
bounds on the spin-independent and spin-dependent 
cross section per nucleon. 



Nonetheless, there are more dimension-6 opera- 
tors that can contribute to the direct detection cross 
section. Namely, 

01 = (X75X)(99), (3) 

02 = (xx)(g75g), (4) 

C3 = (X75x)(975g), (5) 
Oa = (x7M75X)(97^g). (6) 

The operators Oi, ©2, and ©4 are not present 
if parity is a good symmetry of the theory, but 
since parity is badly broken in the Standard Model 
and it could be badly broken in the dark matter 
sector, it is reasonable to include them. If x is 
a Dirac fermion, instead of Majorana, additional 
operators are possible. In particular, there is the 
possibility of a dipole or charge radius coupling to 
dark matter and a vector coupling to quarks Jj l2j . 
Such an operator is quantitatively similar to Oi, 
with the principle difference that it typically couples 
to atomic number Z rather than mass number A. 

These operators in Eqns. (|3])-([6]) are present even 
in the context of the minimal supersymmetric Stan- 
dard Model (MSSM), but there the contributions to 
scattering are typically far subdominant, as they are 
suppressed relative to Osi/sd by additional powers 
of momentum 0(g^/M^) ~ 10^^, or in the case of 
O4 also by velocity suppression u^. Consequently, 
they are usually ignored f3], but see [4]. Moreover, 
even if the dominant operators are zero, because of 
this suppression, these new operators are typically 
negligible in the context of direct detection experi- 
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ments. Thus, even neglecting Osi/sd it niight seem 
unlikely that such interactions would be relevant for 
upcoming direct detection experiments. 

However, this reasoning ignores that these two 
facts often go hand in hand. Osi/sd ^-^e typically 
small when there is a symmetry reason for them to 
be small, in particular, when the DM-nucleon force 
is mediated by a pseudo-Goldstone boson (PGB). 
Because of their shift symmetry, PGBs have 
suppressed interactions. At the same time, PGBs 
are also naturally much lighter than the weak scale. 
Thus, 01-4 are no longer insignificant if the media- 
tor has mass < O(GeV). When we combine this with 
the recent interest in new GeV-scale particles, e.g., 
[5j, and in particular PGBs ^ arising from models 
to explain PAMELA, Fermi, ATIC and HESS, we 
are strongly motivated to consider these scenarios. 

In this paper, we explore a class of dark mat- 
ter models where the scattering is momentum de- 
pendent (MDDM), i.e., where the operators Oi_4 
dominate and are large enough to be observable 
in upcoming direct detection experiments. As we 
shall see, these operators can have a significant 
impact on the spectral shape and the sensitivity of 
various experiments. As an example, we shall see 
that these effects can improve the ability to explain 
the DAMA annual modulation signal while being 
consistent with other direct detection exclusions. 



SIGNALS OF MDDM 

The recoil rates at a direct detection experiment 
can be written as 



dR 



NtitinPx 



2m 



fiv) 



(7) 



where mjsf is the nucleus mass, Nt is the number of 
target nuclei in the detector, p-^ = 0.3 GeV/cm'^ 
is the WIMP density, p is the reduced mass of 
the WIMP-nuclei system, and f{v) is the halo 
velocity distribution function in the lab frame. The 
minimum velocity to scatter with energy En is 
Vmin = \/mAri?i{/2/i2. The rest of the expression 
depends on the scattering's = 2mAr£'^. For SI 
interactions, we have 



SI 



[Zfp + {A~Z)UfF\q'), (8) 



where /p, /„ are respectively the couplings to the 
proton and neutron, and factor is the form factor. 



We take the limit fp = /„. This expression is then 
proportional to the nucleon scattering cross section 
CTp — ^Gpp^pfp, where pp is the reduced mass of the 
WIMP-proton system. For SD interactions, we have 



o'('7^)sD 



2J+1 ^"P^' 



3(9^) + apa^Spn^q^) 



(9) 



where ap,a„ are respectively the couplings to the 
proton and neutron, and the S factors are the 
form factors for SD scattering. The corresponding 
nucleon cross sections are (T(p,„) = ir^F/^(p,n)'^(p,n)- 
The effect of the new operators can be parameter- 
ized simply: 



dER 




dR, 
dEn 



, (10) 



where i indexes the interaction, i.e., SD-proton, SD- 
neutron or SI, and we have included the propagator 
due to a light mediator (j) with mass m^. For the 
benchmark cases, we will take ^ g^, to arrive at 
the simple form 



dR 



MDDM 



dER 




(11) 



We have chosen to normalize the new factors out 
front at a reference value q^^^ = (100 MeV)^, 
a characteristic value for many direct detection 
experiments. For operators 01,02 the exponent 
n = 1, while for O3, n ~ 2. For Oi, the interaction 
is spin-independent on the nucleus side, while for 
the others it is spin-dependent. This form of the 
recoil rate defines the nucleon cross sections (Tp_„ for 
momentum dependent scattering. 04's scattering 
cannot be written in this form, since it has terms 
proportional to the DM velocity. However, we find 
that its spectra is almost identical to standard SI 
scattering, so we neglect it for the rest of the paper. 

MDDM is characterized by a modification of its 
nuclear recoil spectrum. Typically, direct detection 
experiments optimize their searches by going to 
lower energy thresholds, where standard WIMP 
signatures are expected to peak. In contrast, the 
spectrum of MDDM vanishes at zero recoil energy, 
and then can be either peaked or fairly flat over the 
range in question. 

We show in Fig. [l] the spectra of MDDM scenarios 
for the case of SI germanium scattering. As we 
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FIG. 1: Germanium spectra plots with arbitrary normalization versus energy recoil for SI momentum dependent 
scattering of a 100 GeV dark matter mass. Plot a) displays the effect of additional powers of with (f,q^,a,nd 
in solid, long dash and short dash. Plot b) illustrates the effect of on the g* suppressed scenario with = 
(1, .1, .01) GeV in solid, long dash and short dash. 



can see, the spectra differ dramatically from those 
expected for conventional dark matter. The powers 
of suppress the low energy events resulting in 
a peaked spectrum reminiscent of inelastic dark 
matter (iDM) In contrast to inelastic dark 

matter, here the peaking arises without needing a 
coincidence of parameters (specifically, 5 in iDM 
models tuned to the WIMP kinetic energy). The 
spectrum need not be sharply peaked, however, and 
can be broadly spread over a large range of recoil 
energies. The non-trivial propagator of Eqn. ( 10 1 



allows the possibility that events can be suppressed 
for ^ m^, as can be seen in Fig. 1 Finally, 
increasing the dark matter mass shifts the spectra 
to higher energies. Given the possibilities, the lesson 
is that search strategies developed for the simplest 
dark matter candidates are by no means optimal 
for every dark matter candidate. The true dark 
matter candidate may not be one of these simplest 
possibilities, and it is important to cast a wide net. 



EXISTING SEARCHES AND DAMA 

To understand the effects of MDDM, we study 
how these q^ effects can modify the limits arising 
from existing experiments. We show in Fig. |2] the 
limits on interactions mediated by Oi compared 
with limits on standard SI interactions, and in Fig. [3] 
the limits of ©3 when compared with standard SD- 
proton interactions. We follow the procedure laid 
out in Ref. [10] for CDMS HI] and XENONIO [H] 
limits and Ref. for KIMS [T^ limits. Although 
PICASSO [15 and COUPP [15 hmits are compara- 
ble, we only discuss PICASSO limits in what follows. 
Our methods better reproduce their (momentum 



independent) result in the SD-proton case, making 
us more confident that the MDDM limit is realistic. 

Inspecting the exclusion limits of these plots, we 
see important changes with respect to the traditional 
cases. First, consider the SI case with q^ dependence 
(d). While CDMS-Ge and XENONIO remain the 
strongest, KIMS becomes stronger than CDMS-Si 
over much of the parameter space. In the SD-proton 
case, limits from PICASSO are significantly weaker, 
and XENONIO becomes stronger than KIMS in the 
15-25 GeV range. 

These results are easy to understand. Due to 
suppression of low energy events in the MDDM 
scenario, experiments that rely upon low energy 
thresholds (in particular, XENONIO) are weakened 
when compared with others with higher thresholds 
(such as CDMS), which is why CDMS improves 
relative to XENONIO. On the other hand, since 
q^ = 2Afjv-Efl, at a given recoil energy, heavier nuclei 
are preferred by momentum dependent scattering, 
which is why KIMS improves over CDMS Si and 
why PICASSO weakens relative to the other ex- 
periments. Another effect occurs for COUPP and 
PICASSO. In these bubble chamber experiments, 
operation at varying temperature or pressure es- 
sentially integrates the recoil spectrum above some 
threshold. The background from alpha decays is 
known to be a flat spectrum above some specific 
temperature or pressure and is fit to in the data. 
For the broadest MDDM spectra (see Fig. [l]), the 
dark matter signal looks similar to this alpha back- 
ground. Unfortunately, this complicates background 
subtraction and reduces the present sensitivity to 
these models. 

Intriguingly, momentum dependent scattering can 
also modify the interpretation of dark matter ex- 
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FIG. 2: Plots of the SI nucleon cross section ap vs DM mass without (a) and with suppression (b). The colored 
regions show the 68, 90, and 99% CL regions for the best DAMA fit. The 90% exclusions limits are KIMS (orange 
dashed), CDMS Si (red sohd), CDMS Ge (red dotted) and XENONIO (brown dot-dashed). We have taken /p = /„. 



planations of DAMA's annual modulation signal 
and exclusions from other direct detection exper- 
iments. The annual modulation signal |16| I17j . 
originally seen at the DAMA experiment [18], has 
recently been confirmed by DAM A/LIBRA [T^. On 
the other hand, limits from XENON and CDMS 
strongly constrain the simplest dark matter interpre- 
tation of the DAMA experiment: a signal resulting 
from the SI scattering of a WIMP. Explanation of 
the DAMA signal with spin-dependent scatterings 
[201 ET is now also strongly constrained by COUPP 
and PICASSO. 

Fig. [2] shows that adding momentum dependence 
to the SI interactions can only weaken, but not elim- 
inate the limits other experiments put on DAMA 
explanations at the 90% confidence level, at least 
within a Maxwellian halo model. Employing the 
caveats discussed in [10 , alternative statistical tech- 
niques [21] or a non-Maxwellian halo [22] might allow 
a window at low mass when combined with these 
new effects. 

In light of this, we now focus discussion on the 
scenario with the weakest direct detection limits, 
SD-proton scattering, shown in Fig. [Sj These 
plots show that the relative importance of different 
experiments can invert as one adds momentum de- 
pendence, for precisely the reasons described above. 
In fact, the normal SD-proton case [201 [H] which 
is ruled out by PICASSO is allowed for the q'^ 
scenario. Interestingly, these factors are also able 
to improve the fit with DAMA's spectral shape, 
so that there are new masses that can now fit the 



DAMA spectrum. In particular, the mass region 
at ~ 40 — 60 GeV would have normally had a 
shape that was inconsistent with DAMA. Since these 
momentum factors suppress the low energy scat- 
tering, the constraint from DAMA's unmodulated 
event rate [lU] is also weakened, leading to better 
consistency with DAMA's full data set. For these 
plots, we assumed a mediator mass of 1 GeV and 100 
MeV. As we will discuss later in the next section, a 
lighter mediator mass of O(IOO) MeV is more suited 
to generate cross sections of this size. As seen in 
Fig. [sj^c), for this lighter choice of mass, the 10 GeV 
DM mass region survives, but the KIMS limit cuts 
into about half of the higher mass region. We should 
note that our approach to the KIMS limits is not 
aggressive, and does not yield as strong a limit as 
that in |13j. Consequently, a more aggressive limit 
might also be able to exclude this region as well. 
For mediator masses much less than an MeV, the 
momentum independent case is recovered since the 
factors cancel in Eqn. (10), at least in the range 



that the first Born approximation is valid. 

One important point about these scenarios is 
that the expected relationship between direct and 
indirect detection signals breaks down. Typically, 
one assumes that the annihilation proceeds into 
some Standard Model final state. Here, since we 
rely upon the light mediator for our interaction, 
it provides an annihilation channel. Then, if the 
mediator is < GeV in mass, it is natural for it 
to decay dominantly to e.g., electrons, muons and 
pions. The limits from Super-Kamiokande WIMP 
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FIG. 3: Plots of the SD-proton cross section ap vs DM mass without (a) and with suppression (b) and (c), 
where the mediator mass is 1 GeV (b) and 100 MeV (c). The colored regions show the 68, 90, and 99% CL regions 
for the best DAMA fit. The 90% exclusions hmits are PICASSO (gray solid), KIMS (orange dashed), XENONIO 
(brown dot-dashed), and CDMS (red dotted). 



capture are then trivially evaded [531 HI] ■ 



MODEL BUILDING CONSTRAINTS 

For the dark matter scattering to display the 
novel phenomenology discussed here, the scattering 
rate must be dominated by the new operators, and 
not the typical SI or SD coupling. This is not a 
trivial requirement. For comparable coefficients, the 
scattering mediated by operators of Eqns. (|3])-(|6]) 
are suppressed by powers of the velocity or momenta 
relative to Eqns. ([!]) and ([2|. It is possible, however, 
that the coefficients of these new operators are much 
larger than the coefficients of the other operators. 
We discuss this further below. 

If the signal is to be observable at near-future 
experiments, the g^" suppression must be compen- 
sated by a large coefficient for the operator. This 
could be due in part to particularly large couplings 
of a mediator to the Dark Sector or a local over- 
density of the dark matter, but the simplest way to 
get an enhancement is just for the mediator mass 
to be small: dR/dEj^ oc rii^'^. The necessary 
depends on the amount of suppression. If there 
is a single , as in Oi and ©2, a mediator mass 
of a 7710 ~ few GeV, 0(1) couphngs to the DM 
while having Yukawa suppression on quark side, 
one finds a ^ 10^''*' cm^ cross section. This is 
near the interesting region for 02- In the case of 
coherent scattering off of nuclei (as for Oi), this 
would actually already be strongly ruled out by 
CDMS and XENON for a large mass range, a viable 



10^^^ cm^ cross section would require something like 
^ 100 GeV. For suppression, the mass has to 
be O(IOO) MeV to get a ~ 10"'^^ cm-^ cross section. 

In specific cases, large contributions to Oi - O4 
can be expected. If there is a light pseudoscalar 
present that couples both to dark matter and to 
quarks then Oi - O3 can be generated through its 
exchange without generation of either Osi or O so- 
li there is no parity violation, then the expectation 
is that ©3 dominates. On the other hand, if parity 
violation is present, then it is plausible that the 
coefficients of Oi - O3 could all be comparable. 
In this case, it is likely that it would be easiest 
to probe Oi because of its coherent scattering off 
of nuclei. Alternatively, parity violation might be 
confined to couplings in the dark matter sector. In 
this case, pseudoscalar exchange could dominantly 
induce 02- We note that if the light pseudoscalar 
is naively realized as a pseudo-Goldstone, it is 
difficult to sufficiently suppress the contributions to 
Osi- Contributions are induced by exchanging the 
scalar whose vacuum expectation value breaks 
the global symmetry and made the (j) light. 

Interactions with only q^ suppression can con- 
ceivably still dominate over standard interactions 
without significant model-building efforts. The 
simplest example comes from charge-radius or dipole 
couplings to a composite WIMP, whose constituents 
are charged under a new, dark gauge group [1,. 
This generates the phenomenology of Oi straight- 
forwardly, although typically with a coupling to 
instead of A^. If the mediation arises through a 
PGB, Oi_2 can dominate over the scalar exchange 
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as only one vertex will be suppressed by f^, while 
the scalar exchange is suppressed by /|. 

The most challenging model-building comes in 
realizing the q'^ suppressed interaction, without in- 
ducing SI scattering from the accompanying scalar 
mediator. While this seems difficult from the per- 
spective of a standard PGB, it can arise fairly simply 
in SUSY theories. While PGBs are a natural way to 
realize a shift symmetry, such a shift symmetry could 
simply be present in the theory from other origins. 
For instance, in theories with N — 2 SUSY in the 
gauge sector (e.g., [21]), there is a chiral superfield 
partner for every gauge boson. The pseudoscalar 
contained in it possesses a shift symmetry which can 
be thought of as a higher-dimensional gauge symme- 
try, compactified on an S'l /Z2 orbifold. SUSY break- 
ing will make the associated scalar massive. This 
can arise either from i^— term breaking, through 
the operator + ^"f)^ (with X a spurion 

that gets the non-zero i^-term, and (j) ^ superfield 
containing the PGB), or from D-terms, through 
WW^c/) (with W"" the U{1)y supersymmetric field 
strength, and VF^ the supersymmetric field strength 
that gets a non-zero I^-term). Even in the presence 
of this supersymmetry breaking, the pseudoscalar 
remains massless, and will only pick up a mass 
radiatively through diagrams violating both the shift 
symmetry and SUSY. Thus, the scalar contribution 
can be effectively decoupled from the strength of the 
pseudoscalar-mediated interactions. 

Finally, a sizable coefficient for O4 can be gener- 
ated in theories with a light gauge boson that couples 
to the dark matter and mixes with the gauge field 
of the Standard Model. 

There are model-dependent constraints on light 
pseudoscalar mediators. In particular, searches for 
axions can apply. For particles in the GeV range, 
the process T — >■ 7^ is relevant. These branching 
ratios are constrained to be in the range 10~^ — 10~^; 
the precise bound depends on the final state of </> 
decay, (j) — fifl, tt or invisible, [26l - f28] . If the cj) has 
couplings comparable to Standard Model Yukawas, 
the branching ratio is a few x lO^'^ for masses well 
below the T mass. Thus, these bounds constrain the 
(j) coupling to b quarks to be somewhat smaller than 
the Standard Model Yukawa coupling. For lighter 
mediators, depending on the fiavor structure of the 
(f) couplings, K — ^ TTTTi/i may be relevant. The rate for 
the potentially more stringent process — !■ 7r+0 is 
suppressed - a pure pseudoscalar coupling does not 
mediate this process, see for e.g., [29J. The dominant 



contribution to this decay comes from tt — (j) mixing 
[3D], which is model dependent. In cases where the 
pseudoscalar couples only to 3rd generation quarks, 
the Kaon decays are absent; however, the Upsilon 
constraints still apply. Following the procedure in 
|31| , we find that 3rd generation couplings alone can 
generate a detectable rate, as pseudoscalar couplings 
to heavy quarks generate a coupling to GG |32| . 
Incidentally, in general, experimental uncertainties 
(in particular, the light quark contribution to the 
proton spin AE) and parameters like tan /? allow 
a proton dominated coupling to be generated. In 
the minimal case of 3rd generation couplings, the (j) 
decays to two photons with a decay length cr ^ 1 m. 
While couplings to leptons are not necessary to 
implement the scenario at hand, if the mediator 
couples with Yukawa strength to the muon, requiring 
that the magnitude of the contribution to the muon 
g — 2 is no larger than the current discrepancy 
between theory and experiment |(5a^| < 290 x 10~^^ 
enforces > 300 MeV |?9j . Finally with couplings 
to electrons, it is possible to search for e+e^ — > (pj 
for either invisible or electron decays of (j) [55] . 
However, for pseudoscalar (/>, suppression by the 
electron yukawa coupling makes the production rate 
below the projected sensitivities |33j . 



CONCLUSIONS 

As the sensitivity of new dark matter direct 
detection experiments continues to increase at a 
rapid pace, the ability to test for new scenarios 
for dark matter will grow simultaneously. Present 
experiments are optimized to search for WIMPs with 
signals that peak at low nuclear recoil energies. In 
contrast, models with momentum-suppressed inter- 
actions (MDDM), have spectra that peak at inter- 
mediate energies, thus changing the expected signals 
and relative strengths of various direct detection 
experiments. While interesting scenarios, specifi- 
cally inelastic dark matter, have been proposed with 
spectra that peak at high recoil energy, we find that 
the scenarios with this feature are more ubiquitous 
than previously thought. In models with new light 
vectors or pseudoscalars, momentum dependent in- 
teractions can be large, and this phenomenology can 
be present. 

A simple parameterization captures much of the 
relevant phenomenology. Specifically, one can re- 
place dRi/dEji with (q^qq)"^ dRi/dEu, where i in- 



7 



dexes the interaction type and gioo is the momentum 
transfer in units of 100 MeV. While additional 
features can arise at low mediator masses, this pa- 
rameterization is sufhcient to reproduce the peaking 
in the spectrum, and provides a convenient way 
to compare different experiments. In analyzing 
the presently allowed parameter space in this way, 
we find that momentum dependent couplings can 
open allowed ranges of parameters for DAMA with 
dominantly spin-dependent proton couplings, if ac- 
companied by an additional suppression. 

Whatever model of dark matter nature has chosen 
to realize, it is important to be cognizant of the wide 
range of possible phenomenology, so that possible 
signals are not missed or attributed to backgrounds. 
The framework of MDDM provides motivation, and 
a prescription to study and constrain these models 
in the future. 

Note added: As this paper was being finished, we 
became aware of |34j . which appeared in the arXiv 
and discusses momentum dependent interactions 
arising from the couplings to one or more new gauge 
bosons, and their ability to explain DAMA from 
spin-independent interactions. 
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